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We discuss cosmological models involving homogeneous and isotropic Yang-Mills (YM) fields. Such 
models were proposed recently as an alternative to scalar models of cosmic acceleration. There exists 
a unique SU (2) YM configuration (generalizable to larger gauge groups) whose energy- momentum 
tensor is homogeneous and isotropic in space. It is parameterized by a single scalar field with a 
quatric potential. In the case of the closed universe the coupled YM - doublet Higgs system admits 
homogeneous and isotropic configurations too. While pure Einstein- Yang-Mills (EYM) cosmology 
with the standard conformally invariant YM action gives rise to the hot universe, Einstein- Yang- 
Mills-Higgs (EYMH) cosmology has a variety of regimes which include inflationary stages, bounces, 
and exhibits global cycling behavior reminiscent of the Multiverse developed in time. We also 
discuss other mechanisms of conformal symmetry breaking such as string-inspired Born-Infeld (BI) 
O ' modification of the YM action or field-theoretical quantum corrections. 
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Q 1 I. INTRODUCTION 

CT 

5— i [ With discovery of inflation as solution of the horizon and flatness problems in cosmology [1] it became widely 
accepted that, apart from gravity, some other homogeneous and isotropic fields have to be present at cosmological 
scale which mimic variable cosmological constant. Traditionally this role is attributed to a single scalar field, the 
inflaton, or several scalar fields Q|. Modern theories provide various candidates for inflaton varying from Higgs field 
of the standard model to more hypothetical moduli fields originating from compactified supergravity/string theory. 
Still, physical nature of the inflaton is far from being uniquely understood, and the choice of the inflaton potential 
remains mostly phenomenological. Similarly, popular current models of dark energy Q involve scalar fields with 
rather exotic properties (quintessence, K-essence, phantom) whose physical nature is far from clear. Moreover, no 
massless elementary scalar field was observed experimentally so far. Meanwhile, vector fields certainly do exist (and 
are masseless before the spontaneous symmetry breaking) being basic ingredients of the Standard Model and its 
generalization. Therefore an idea to use vector fields instead or together with scalar ones to model inflation and dark 
energy seems to be appropriate. In fact, the model of inflation driven by vector field was suggested long ago [H, but 
remained unnoticed by cosmologists until recently when it was revived in the context of the dark energy problem [f| . 
Formation of YM condensates in superdense matter was discussed long ago by Linde [6|. 

There are two major reasons why vector fields were not welcome in cosmology, apart from their relative complexity. 
Spatially homogeneous configuration of a single (Abelian) vector field evidently can not be isotropic. Therefore, 
in order to fit the Friedmann-Robertson- Walker (FRW) cosmology one has to introduce (at least) a triplet of vector 
fields ensuring isotopy of the total energy-momentum tensor. Another problem is conformal invariance of the standard 
classical YM lagrangian, which leads within the FRW cosmology to the equation of state w = P/e = 1/3, identical 
with that of the photon gas. Therefore, the solution for the coupled EYM system will be the hot Universe driven by the 
cold classical matter field Meanwhile, for an accelerated expansion one needs the equation of state —1 < w < —1/3, 
so the conformal invariance must be broken. 

However, the first problem is avoided in the non- Abelian case: the SU(2) triplet of YM fields admits an essentially 
non- Abelian configuration (with non-zero commutator of the matrix-valued potentials) whose stress-tensor exhibits 
three-dimensional homogeneity and isotropy. The second problem (conformal symmetry breaking) can be overcome 
passing to effective actions which account for quantum corrections either in the context of gauge theories or string 
theory. Various attempts to use YM vector fields in constructing dark energy models thus were undertaken during 
recent years Q. 
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II. SU(2)-DRIVEN FRW COSMOLOGY 

Consider the FRW interval in the conformal gauge 

ds 2 ^a 2 ^)(dr, 2 -dll) (1) 

where k = — 1, 0, 1 for open, fiat and closed spatial geometry. As it was shown for k = 1, by Cervero and Jacobs 
Henneaux [l(| and Hosotani and for all k by Gal'tsov and Volkov [tJ, the following configuration 

E? = f5?, B<; = {k-f 2 )5t, (2) 

parametrized by the single scalar function f(rj) of the conformal time gives rise to homogeneous and isotropic energy- 
momentum tensor. The YM lagrangian density in the conformal frame then reads 

so the electric part corresponds to the kinetic term, while the magnetic part — to the potential term in the action. 
The effective scalar field / is dimensionless (we use the units H = c = 1), so in spite of similarity with the scalar 
field theory with the quatric potential, dependence of the lagrangian density on the scale factor a is different: the 
energy scales as 1/a 4 which is characteristic for the conformal field. The energy- momentum tensor is traceless and 
the equation of state is 

P = e/3 with e =-i 7 [/a + (A ; -/ 2 ) a j. (4) 

Thus one obtains the hot Universe driven by non-thermal matter Q: our classical YM configuration perfectly mimics 
the photon gas. 

The YM equations reduce to equations of motion of a fictitious particle in the potential well 

W k = (k- ff (5) 

which is especially interesting in the closed case k = 1 when it is the double well potential (Fig.l). Two minima 
correspond to neighboring topologically different vacua. We therefore observe that gravity lowers the potential barrier 
between topological sectors to a finite value, similarly to the Higgs field. Physically this similarity is due to attractive 
nature of both (contrary to repulsive nature of YM). If the energy is less than the height of the potential barrier, 
the particle oscillates around a single vacuum, when it is above the barrier, oscillations between different vacua are 
unsuppressed. 

In the flat and open cases the potential has one minimum at / = with W = in the flat case k = and W — 1 
in the open case k = —I. so there are no topological effects. 
Computing the Chern-Simons 3-form 



uj 3 = — Tr \Af\dA- —A A A A A] , (6) 



satisfying the equation 



du 3 = — jTrF A F, (7) 

one finds that it is non-trivial in the closed case k = 1, giving the winding number of the map SU(2) — > S 3 : 

N C s= [ u; 3 = i(/ + l) 2 (2-/). (8) 
Js 3 4 

The vacuum / = — 1 is topologically trivial: Ncs — 0, while the vacuum / = I is the next non-trivial one with 
Ars 1. 

Generalization of the above ansatz for SU(n) and SO(n) gauge groups and further classical and quantum properties 
of EYM cosmological solutions were considered in a number of papers (l2l . [l3l in the 90-ies. Behavior of small 
perturbations in cosmologies with vector fields was discussed more recently in |14j . 



Closed FRW Universe 
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A. Cosmological sphaleron 

In the closed case there exists a particularly simple configuration / = (/-particle sitting at the top of the 
barrier between two vacua) which in analogy with the sphaleron in the Weinberg-Salam (WS) theory was called the 
"cosmological sphaleron" [15j . It is worth noting that the localized particle-like EYM solutions similar to the WS 
sphalerons exist as well which are asymptotically flat regular particle-like solutions of the EYM equations discovered 
by Bartnik and McKinnon (for a review and further references see [lH). The repulsive YM stresses in these objects are 
compensated by gravity instead of the Higgs field in the WS sphalerons. Creation and decay of sphalerons generates a 
transition of the YM field between topological sectors, and it is accompanied by the fermion number non-conservation 
in presence of fermions [16] . The cosmological sphaleron has the topological charge Ncs = 1/2 like the sphaleron in 
the Weinberg-Salam theory. 

The equation of motion of the /-particle 

/ = 2/(1 - f) (9) 
is solved indeed by / = 0, this correspond to the total energy 

f 2 + (i - Z 2 ) 2 = i. (io) 

The YM field in this case is purely magnetic. A more general solution with the same energy describes rolling down of 
the /-particle (the sphaleron decay) |15j : 

- * an 



cos 



Rolling down to the vacuum w — 1 takes an infinite time, while the corresponding full cosmological evolution is given 
by 



UirG , . 

— 2-sinr? (12) 



and takes a finite time. Thus the cosmological sphaleron is quasi-stable. This conclusion is not modified if a positive 
cosmological constant is added. 



B. Instantons and wormholes 



An homogeneous and isotropic EYM system has interesting features also in the k = 1 space of Euclidean signature 
which is invoked in the path- integral formulation of quantum gravity. Actually, when |/| < 1, transitions between 
two topological sectors can be effected via underbarrier tunneling described by instanton and wormhole Euclidean 
solutions. In the Euclidean regime the first integral of the equations of motion reads: 



f 2 - (f - 1) 2 = -c, 



(13) 
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where C is an integration constant. Instanton corresponds to C = 0, it describes tunneling between the vacua / = ±1. 
It is a self-dual Euclidean YM configuration for which the stress tensor is zero, therefore a conformally flat gravitational 
field can be added just as a background. 

Tunneling solutions at higher excitation levels < C < 1 are not self-dual. In flat space-time they are known 
as a meron (C = 1, the Euclidean counterpart of the cosmological sphaleron, Nqs = 1/2) and nested merons 
< C < 1, 1/2 < Ncs < 1. The energy- momentum tensor of the meron is non-zero, and in flat space this solution 
is singular. When gravity is added, the singularity at the location of a meron expands to a wormhole throat, and 
consequently, the Euclidean topology of the space-time transforms to that of a wormhole. Topological charge of the 
meron wormholes is zero, the charge of the meron being swallowed by the wormhole [13] ■ The total action of these 
wormholes diverges because of slow fall-off of the meron field at infinity, so the amplitude of creation of baby universes 
associated with the Euclidean wormholes is zero. However, when a positive cosmological constant is added (inflation) 
the action becomes finite due to compactness of the space. Such solutions can be interpreted as describing tunneling 
between the de Sitter space and the hot FRW universe. 

Adding the positive cosmological constant, we will get similar first integrals both for the /-particle and the 
cosmological radius: 

f 2 - (f - I) 2 = -C a 2 + (Aa 4 /3 - a 2 ) = -C/(e 2 m 2 Pl ). (14) 

Solutions describe independent tunneling of / and a with different periods Tf, T a depending on the excitation level 
C. To be wormholes, they must obey a quantization condition njTf — n a T a (l8| with two integers. However, for a 
specific value of the cosmological constant 

3 

A = -m 2 Pl , (15) 

it was found (l9j that T a = T f for all C e [0, 1]. In particular, for C = 1 (the meron limit) the radius a becomes 
constant (Euclidean static Einstein Universe). For C / 1,0 the solutions describe creation of the baby universes 
(which was invoked in the Coleman's idea of the "Big fix"). Remarkably, under the above conditions, the total action 
(gravitational plus YM) is precisely zero [l9j: 

Sym + S gr = 0. (16) 
Thus, the pinching off of baby universes occurs with unit probability. For later work on EYM wormholes see [2oj | . 

III. EINSTEIN- YANG-MILLS-HIGGS COSMOLOGY 

Consider now the EYMH action with complex doublet Higgs: 

S = J {-^R-lF2 v F^ + l(D^D^-±(&$-v 2 ) 2 }v^d 4 x, (17) 

where 

= + eA^Ta®. (18) 
It is easy to see that in the case of spatially closed FRW cosmology an ansatz for Higgs 

$ = h(t)e im U<S>o, U = e 2xT °"°, $^$ = 1, (19) 
is compatible with homogeneity and isotropy of the full EYMH system. Indeed, with the YM parametrization 

eA^T a =^—Aud tl U- 1 , (20) 
the covariant derivative reduces to the ordinary one 

= ^rA<F (21) 

Thus an homogeneous and isotropic cosmology does exist for Higgs in the fundamental representation. On the contrary, 
the triplet Higgs leads to anisotropic cosmology [2l| . 
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The Higgs phase rotation factor e'^*), which is compatible with the desired symmetry, in the fiat space leads to 
infinite energy and must be omitted, but for the closed FRW cosmology its contribution is finite, so we may keep it. 

The EYMH action contains three different mass parameters: the Planck mass, the mass of the W-boson, and the 
Higgs mass 



M H = V\v. 



We then rescale the Higgs function h —> hMpi and introduce the dimensionless parameters 



Mpi = — =, M w = ev, 
v G 



Ah 



w 



eMpi Mpi 



M 



1R 



M w e 



(22) 



(23) 



Finally, keeping in mind that the pure EYM system has a natural lengths scale I = l/(eMp\), we present the metric 
in terms of the dimensionless lapse and scale functions N, a: 



ds 2 = I 2 {-N 2 dt 2 + a 2 [d X 2 + sin X 2 (d8 2 + sin 2 9dip 2 )] } 



(24) 



Then substituting the ansatze (|19l20j) and integrating over the three-sphere we obtain the one-dimensional reduced 
action 



Si 



3 

8^ 



aN - 



N 



3f 2 a K 2 a A h 2 Z 2 a 3 



2N 



2N 



2N 



Vf - V h - Mnt 

3 



dt 



Vf = mP 2a ^ , V h = ^{h 2 - a 2 ) 2 Na\ V int = -^{f + l) 2 Na 

where we have omitted an overall factor 7r/4 and the total derivative in the scalar curvature term. 
The equation for the Higgs phase rotation variable £ 



(25) 
(26) 



d_ 

dt 



h 2 a^ 
N 



= 



(27) 



can be integrated 



V2jN 
h 2 a 3 ' 



(28) 



where j is the integration constant. Then it is easy to check that the remaining dynamics can be derived from the 
action in which the ^-kinetic term is replaced by the potential term 



V, = -J—. 

3 Nf 2 a 3 



Variation with respect to N leads to the constraint equation 



a{a 2 + 1) = 



87T 



(29) 



(30) 



and it is convenient to fix the gauge N = 1 afterwards. For the energy density we then get 

e = T f + T h + V f + V h + V int + Vj , 

where the first two terms are kinetic: 



(31) 



2„3 



-h a 



T 



f 



r 



(32) 



and the remaining are four potentials: 

3(/ 2 - l) 2 



Vf 



2a 



v h -^n. 



-(h>-a 2 ) 2 a\ Vint = h 2 (f + l) 2 a, \ , j^. 



(33) 
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The acceleration equation is obtained by variation of the action with respect to a with account for the constraint: 



2- 47T 

a a = — — (£ + op), 



(34) 



leading to the following expression for pressure: 



rji T/ . T r . Tf Vf Vint 

P = T h - V h + Vj + -± + -± - —, 



(35) 



and therefore, 



a 2 d = -^-(2T h + 2Vj + Tf - V h + V f ). 



(36) 



The field equations for the YM and Higgs scalar functions read 



f + a -f 

a 



-^ (/ + l)-^(/ 2 -l)/. 



(37) 



From the acceleration equation is clear that kinetic terms always produce deceleration, while the potential terms are 
partly accelerating. In the transient regimes when different potential terms dominate, one has the following equations 
of state: 



Dominant potential 


w = p/s 


type 


Higgs V h 

Higgs phase rotation Vj 
YM potential V f 
Interaction term Vint 


w = —1 
w = 1 
w = 1/3 
w = -1/3 


cosmological constant 
stiff matter 
radiation 
string gas 



A. Standard model scale 



Behavior of solutions essentially depends on the parameters a, (3. Consider first the scales relevant to the Standard 
model. In this case a ~ 10~ 17 and f3 is of the order of unity. The corresponding dynamical scale if far from the 
Planck scale and matter contribution to the action is of the order of a 4 /? 2 . The corresponding values of the scale 
factor must be of the order l/a 2 (3. Then the main contribution to evolution of the scale factor will come from the 
Higgs potential Vh- Typical regimes of the evolution of the Higgs scalar correspond to motion near the extremal 
points of the potential: the minima \h\ = a and the local maximum h = 0. In the first case one observes small 
oscillations, in the second - slow rolling down from the top of the potential barrier. While the frequency of oscillations 
is proportional to (3, the rolling down velocity is much less, namely of the order of a 2 (3. A substantial variation of the 
scale factor will correspond to time intervals of the order of l/a 2 /3, so the cosmic acceleration will be proportional to 
the average in time values of potential and kinetic energy of the scalar field < Vh — 2Th >■ For harmonic oscillations 
< Th >~< Vh >, so the acceleration will be negative. The rolling down regime is exponential, so the characteristic 
time T ro u will depend on initial conditions as — ln/i(0). During the rolling time the scale factor will exponentially 
grow up with the Hubble parameter proportional to a 2 (3. Thus an exponential expansion will be insignificant unless 
we choose — ln/i(0) > l/a 2 f3. 

It is known that in inflationary models with power-law potential the oscillating universe regime is possible when 
during contraction phase the amplitude of scalar field oscillations grows up and as a result the scalar field climbs close 
to the potential top. Then the rolling down regime follows with the corresponding expansion of the scale factor. When 
rolling down terminates, the oscillations starts again and one enters a new cycle. But with the YM field present, the 
probability to hit the relatively small region of the phase space for triggering the rolling-down regime decreases, since 
the cumulating chaotic oscillations and interaction with the YM component deviate substantially the phase trajectory 
from reproducing precisely the previous cycle. Therefore for small a the cycling behavior demands fine tuning of the 
initial data for realization of the cycling universe. Otherwise, the evolution will terminate with a collapse to a point. 

As about the YM component, its amplitude will oscillate with a period inversly proportional to the average value 
of the scalar field <h> (i.e. of the order a -1 ) around the minimum / = —1, since oscillations will be governed by the 
interaction potential V m t, while the YM potential Vf will be negligibly small. In the case when the scalar field sits 
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at the top of the potential barrier (what is possible if Vj — 0), the YM filed will oscillate in the potential Vf. The 
oscillation period will then be proportional to the value of the scale factor, but for sufficiently fast increase of a the 
frequency becomes imaginary, describing an exponential relaxation of the YM field to the vacuum value. The scale 
factor itself in this special case will be either exponentially growing, or collapsing depending on initial conditions, 
This can be shown analytically. Indeed, for h — the non-zero potential energy of the Higgs field will play the role 
of the cosmological constant equal to A = a 4 /3 2 /4. As a result, independently of the dynamics of the YM field the 
equation of state will be 



Substituting this into the continuity equation 



after trivial integration we obtain 



£ - 4Aa 3 

P= o • ( 38 ) 



i + 3-p = 0, (39) 
a 



C 2 

e = Aa 3 + — , (40) 
a 



where C - is the integration constant. For instance, for the sphaleron configuration / = h = with no field dynamics 
the energy density is given by 

e = Vh + V, = ±£ + l, (41) 

which results in C 2 = 3/2. Using the value (|40p for the energy density, we can solve the constraint equation for the 
scale factor: 

/2-k [2tt 

d 2 + 1 = 7 2 a 2 /4 + 4?7 2 /a 2 , where 7 = 2a 2 /3Wy, r]=CJ—. (42) 

A simple solution arises when the right hand side of this equation is a full square. For this to happen, one has to 
impose the following relation between the parameteres: 

777 = 1. (43) 

Taking the square root of the constraint equation, one obtains 

= 70/2-2/(70). (44) 

Then the solution for the scale factor will read 

/ 1 / 1 \ 

(45) 

r,n,e a(0) = ao- Note that the constraint equation is invariant under the time reflection, so if a(t) is a solution, then 
a(— t) will also be a solution. So for brevity we do not write ± in the Eq. ([4"4"| as well as in the solutions for a(t). 

Depending on the values of the parameters, one observes either cosmological singularity, or inflation. For ao = 2/7 
the solution will be the static universe a(t) = a : in this case the negative pressure of the scalar field is exactly 
compensated by the positive contribution of the YM field. 

If the right hand side of (l4"2"j) is not a full square (which can be expected in our Standard model scales case since the 
relations (14"3"|) does not hold in view of a <C 1), then an analytic formula for the time dependence of the scale factor 
is more complicated: 




_/ 1 _ 7 2 ?? 2 + 

7 4 C 7 2 



a = \l— (!-7V) + — + CeT. (46) 



Here the integration constant was introduced as C and not as a(0). Evidently, for small a the solution will be a 
growing exponent with the power coefficient proportional to a 2 . Coming back to dimensioned time parameter will 
get an exponential of a 2 /Tpi ~ I0 9 sec _1 . Comparing this with an observed Hubble constant 10~ 17 sec _1 we get the 
needed value of the mass of the scalar field 10 _30 Afp; = 10 -2 eV/c 2 . 



B. Planck scale 



Now consider the case of a and /3 of the order of unity, so all the quantities are of the Planck scale. While in the 
case a 2 (3 <C 1 the behavior of the system was determined by the fixed points of the Higgs potential, now we have to 
find the fixed points of full system of equations (157)) and (j3"u) . 

Equating to zero the time derivatives of the variables we get the following system : 

~ 3h ;(f + l) 2 + hf(h 2 -a 2 ) = 0, 



2a 2 

¥(/ + i) + % 2 -i) = o, 



^-{^-c?f-^{f-lf - °- ( 47 ) 

Denoting the variables collectively Q — {h, f, a} we will get a solution depending on parameters Q — Qo(a, (3). This 
solution has to satisfy the constraint equation ([50)) : 

a (a,P)=^-e(a,P). (48) 

Not every solution of (|47[) does it: e.g. the vacuum state h — ±a, / = — 1 does not. As a result, the physical fixed 
points of the system of equations are realised in the parameter space only on the curve given by the Eq. (|48)) . 

Obviously the solution / = h = satisfies the first two equations of the system (|47[) . while from the third equation 
we find the complete solution: 



h = 0, f = 0, aa = 2/y/3. (49) 

He we denoted /3 = /3^/8/3 to simplify further relations. This static solution was already found in the previous section 
as ao = 2/7. 

Finally, the last non-trivial solution of the system (|4"Tj) reads: 

,2 , 0-1 

an = j = a— , aa = — . (50) 

/3 + 1 y/3 + i 



In the limit (3^-1 + this solution coinsides with the previous one. For j3 < 1 the solution of the type (|50|) does not 
exist. One can notice that the parameter a us just the scale one. 

Let us try to give qualitative description of the fixed point corresponding to the solutions described above. For the 
first solution we have the following configuration. The field h sits at the local maximum of its potential Vh' ho = 0. 
This switches off the interaction between the scalar and the YM components: I^nt = 0. Correspondingly, the fixed 
points of the YM field will be extrema of the potential Vf. These are ±1, 0, from which only /o = satisfies the 
constraint equation. 

In the second case the fixed point for h will be not the maximum, but the minimum of the potential. This is not the 
minimum h = ±a, since interaction with the YM field shifts the point of minimum in such a way that < ho < a. 
Similarly, for the YM field the extremal points 0, 1 turn out to be shifted: the local maximum to the right, and the 
minimum to the left according to the relation f max /min = (1 i \/T— a 2 h 2 )/2. From the constraint equation we get 
that for / < 3, the fixed point is the local maximum; for I — 3 the maximum and the minimum coincide forming the 
inflection point; with further increasing I the fixed point will be the local minimum. 

Finally, the constraint equation shows in which region of the phase space the system can be depending on parameters, 
or, conversely, which should be the energy of the scalar field for the desired regime of evolution associated with the 
motion near one or another fixed points. Therefore for two families of solutions (|49[) and (|50[) . respectively, we find 
the following relations: 

* = - p=, a 2 = w =. (51) 

To describe evolution of the system in the vicinity of fixed points one has to linearize the equations of motion 
around them. Excluding the momentum variables, we obtain the matrix equation: 

SQ = MSQ, (52) 
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where the matrix M is obtained by differentiation of the equations of motion and the Friedmann equations over 
Q = {h, f, a} and substituting the values of the variables in the fixed point. At the same time we have to take into 
account the dependence between the parameters imposed by the constraint equation. Then we obtain for M(l); 



while in the second case: 



M 9 



Mi 



3Q3 2 -1) 
32tt 



f(/3-l) 




3V2(/3+l) 5/2 3V2(0+lf yi^T \ 
128tt 3 / 2 ^ 2 512tt 2 /3 2 

3(/3-3)(j3+l) 3(ff+l) 3 ^ 2 (/-l) 

167T/3 ~~~ 32tT 3 / 2 ^ 2 



(53) 



V 



13-1 
2 



Q3 + l) 3 / 2 (ff-l) 
8v^/3 2 



16TT/3 2 



(54) 



It is easy to explore the case j3 < 1, when only the matrix Mi is relevant. For the scale factor the fixed point will 
be the focus by h and nodes by / and a describing an exponential expansion or contraction of the universe. For this 
we obtain the following critical value of the scale factor ao — \/8tt. when the inflationary potential Vh is compensated 
by the deflationary potential Vf in the equation. 

In order to determine the system behavior in the general case, one has to find the eigenvalues of the linear system. 



det(M 1)2 - fi 2 I) = 0. 



(55) 



where the eigenvalues fx of the full system in the six-dimensional phase space enter squared into the three-dimensional 
system which is left after eliminating momenta. 

The matrix Mi is diagonal, so that the eigenvalues are obvious. The eigenvalue for the scalar field is/i^ = 3/4(/3 — 1). 
For P > 1 we have a node. For /3 < 1, the eigenvalues become pure imaginary, that is the node is transformed to 
focus. This means that the quadratic in h interaction potential Vi nt dominates over the quadratic term of the Higgs 
potential, so the effective potential at the point h = has a minimum, but not a local maximum. Eigenvalues of the 
gauge field and the scale factor are obviously always real. On the phase portrait this fixed point is a node, and one 
therefore has expansion or contraction of the universe. Correspondingly, the system is unstable in the vicinity of fixed 
point parameterized by the equations (|49|). 

Eigenvalue of the matrix M 2 can also be found analytically since the corresponding equation (|55|) is cubic in terms 
of /i 2 . Their explicit form, however, is too long, so we do not present it here. One root is always real, as it should 
be for the cubic equation. It is negative and the corresponding fixed point is focus. Two other roots are complex 
on the interval /3 £ [1.8, 7.5]. Their real parts are equal, while imaginary parts differ in sign describing different 
orientation (left/right handed) of phase trajectories. For 1 < f3 < 1.8 the roots are real and positive, while for j3 > 7.5 
- negative. In the first case one has nodes, in the second — elliptic orbits. Thus, on the interval /3 > 7.5 the fixed 
points parameterized by (|50l) . are focal points in the full six-dimensional phase space, so the system is stable in the 
vicinity of these points. 

From this analysis it follows that one of the main features of the behavior of the system on Planck scale is that 
apart from unstable solution of the type of static universe, when small deviation trigger long cycles of expansion 
or contraction, there are quasi-static regimes characterized by dynamical equilibrium between the scale factor and 
matter, which are stable against small perturbations. 

Finally, we discuss the role of the Higgs phase rotation parameter j. When j ^ 0, one observes the regimes typical 
for the second family of fixed points (while the first family can not be realized since now ft ^ 0). This additional 
parameter can be used to set the system into the state close to stationary points. Without it, we have some prescribed 
dependence between the other parameters, a(/3), for every fixed point, which may not be satisfied for the particular 
theory chosen. With j ^ 0, one has more flexibility in the parameter space to impose the desired evolution regime of 
the system. 



C. Numerical experiments 



Numerical solutions confirms the above qualitative considerations. We are particularly interested in illustrating the 
system behavior in different regimes described in the previous subsection. On Figs. (jTH3) we present the solution 
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exhibiting the existence of quasi-stable state. The system enters this state with negative and leaves it with positive 
cosmic acceleration, with concrete values depending on initial conditions. The Fig. §3§ illustrates transition from the 
quasi-stable local minimum near zero to a true minimum h = a, after inflationary expansion of the scale factor. These 
figures show behavior of the system with Planck scale parameters, when interaction with the YM field is sufficiently 
strong to change substantially vacuum state of the scalar field. 

The cyclic evolution is presented on Fig. (j4]). For this solution we have chosen the parameter values such that 
the periods of the expansion cycles are comparable to the period of the oscillations of the scalar field. In this case 
interaction between scalar and YM fields is not strong, but still the YM field plays crucial role at the moments of 
bounces. Since the oscillation frequencies of h and / are different, the vaule of the YM field at the next bounce may 
differ significantly for that at the previous bounce. With high probability the subsequent cycles differ from each other 
with the YM field values being chaotically distributed at the moments of bounces. This is the main difference of the 
large time scale behavior of the EYMH system as compared with the pure scalar cosmology. 

Finally, on Fig. ([5]) we illustrate the dynamics of the EYMH universe for asymptotically small value a = 10~ 6 . 
When the oscillation amplitude of the scalar field is relatively big (for this figure h(0) = 10~ 2 ), the potential term 
exceeds the kinetic term leading to weakly accelerating regime. But with growing scale factor the scalar field amplitude 
decreases and expansion converts to contraction. 

D. EYMH hybrid inflation 

As was described in previous subsection, the EYMH system demonstrates a wide variety of behavior in different 
regimes. But the most interesting among them seems to be the next one. The EYMH system may be considered as 
a very natural variant of the hybrid inflation, introduced by Linde [22l |. The main idea of the hybrid inflation was to 
add a massive scalar field / to the Higgs field h so that the total potential reads 

V(f,h) = ^-f + ^fh" + ^-v^. (56) 

In this case inflation consists of two stages. When / is larger then the critical value f c — Xv 2 /X', the Higgs field is 
trapped on the top of its potential due to the interaction between the scalar fields. During this stage the field / slowly 
rolls under f c and then the second stage starts — the roll of h. Inflation ends when h reaches its true minimum |t>|. 

In our system there is a similar trapping potential V; n t . Mention that it does not enter the Friedmann equation for 
acceleration, therefore the inflationary potential will be only Vh- Next, Vint depends on the gauge field and the scale 
factor. The gauge field changes slowly (especially if one sets it in the minimum / = 1 of its potential Vf), but the 
scale factor in the denominator grows exponentially, therefore the first stage when h is trapped will be rather short. 
But during the first stage the scalar field will be lifted closer to the top of the potential Vh, what implies a prolonged 
second stage of inflation with the slow roll of the Higgs field. In other words, large inflation occurs when the initial 
value of the Higgs field, ho, is strictly chosen to be in vicinity of the top: ho ~ 0. But this condition can be weakened 
significantly due to the trapping potential which will automatically prepare the Higgs field in the right position before 
the rolling down. 

As was mentioned above, the gauge field plays crucial role in the trapping of the Higgs field: it switches the trapping 
on, when fo = 1, and off, when fo = — 1. In the last case the Higgs field will fall freely, reproducing the standard 'slow 
roll' scenario. The initial value of the scale factor ao affects the duration of the trapping: when a reaches the critical 
value a c ~ \/6/(a0) (what happens rather fast due to the exponential growth of a), the trapping will end. 

Now let us turn to the numerical calculation of the number N of e-folds during this hybrid inflation. On the Fig. © 
one can see the dependence of the N on the initial value ho of the Higgs field. There are three plots: when trapping 
is switched off by the gauge field fo = — 1; N+, which was calculated for the trapping with fo — 1 and ao/a c = 1/5; 
their ratio N + /N_. The valued N± are normalized by the number N = 60 of e-folds during the real inflation of our 
universe (23| . One can see that even for rather weak trapping with ao being just five times smaller than the critical 
value, there is a 20 — 40% gain of e-folds. This gain grows with the increase of ho, as expected, because the lift of h 
to the top of Vh is significant when ho is far from the top. Also the needed number of e-folds N ~ 60 can be obtained 
in twice wider range of ho when the trapping is on. 

So switching on of the trapping due to the interaction with the gauge field may really enhance inflation. On the 
Fig. we can see how this amplification will depend on the ratio In (ao/a c ). For chosen parameters the absence of 
trapping when ao S> a c gives us about thirteen e-folds. But the choice of several e-folds smaller ao (so that ao <C a c ) 
will allow us to gain a large total number of e-folds. Actually there is even a divergence in the e-folds number when 
the trapping potential drives the scalar field exactly on the top of the potential Vh with infinite inflation. Of cause 
this trapping scenario is significant when j = 0, and the number of e-folds greatly decreases with the growth of j. 
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We conclude this section with the following remarks. Coupled Yang-Mills-Higgs dynamics for closed FRW universe 
gives rise to new interesting evolution types which include transient regimes of cosmic acceleration, bounces and cyclic 
evolution. Presence of the YM component and Higgs phase rotation parameter changes substantially the dynamics 
of the universe at small scale factors. The interaction with the gauge field holds the Higgs field near the top of the 
potential, and the balance of accelerating and decelerating potentials of scalar and vector fields can freeze the scale 
factor. The phase rotation acts in reverse, making the dynamic of the system to be more similar to the scalar field 
with power-low potential, but its kinetic nature leads to the opposite sign in the potential at very small distances. 
Also it increases the number of free parameters, which can be useful in quantitative analysis. 

One intriguing feature is possibility of an infinite sequence of cycles whose parameters change chaotically due to 
evolution of the YM component. This resembles the Multiverse models [lij realized as sequence of universes with 
different parameters in the ultralarge time scale. 

The system of interacting Higgs and YM fields can be considered as a very natural candidate for the hybrid inflation 
scenario, which is richer then a standard slow roll inflation with a scalar field. The particular feature of the model 
is that due to the vector nature of the gauge field, whose energy density depends on the scale factor, the EYMH 
inflation also inherits this dependence. There can be a large inflation in a small universe, and a small inflation in a 
large universe. This looks quite similar to the current views on the evolution of the universe with large initial inflation 
and slow late-time acceleration. 



IV. NON-ABELIAN BORN-INFELD (NBI) 

Open string theory suggests the following generalization of the Maxwell Lagrangian (applicable to any dimensions) : 



?2 



L = ^ LP det(g^ + F~/fi) - , (57) 

j3 being the critical BI field strength (j3 — l/2%a' in string theory). In four dimensions this is equivalent to 



In the non-Abelian case the strength tensor is matrix valued, and the prescription (Tseytlin [25]) is more 
complicated: the symmetrized trace, which is calculated expanding the Lagrangian in powers of F£ v T a (T a being 
the gauge group generators), then symmetrizing all products of T a involved and only afterwards taking the trace. 
Symbolically this is given by the expression 

L str = ^Str (y~ det( ffMV +F MV //3) - y/=g\ , (59) 



47T 

but actually this is a useful form if one is able to perform a subsequent resummation. Fortunately, this is possible in 
the closed from for the homogeneous and isotropic SU(2) YM field and the metric ds 2 — N 2 dt 2 — a 2 dl 2 leading to [2^ : 

at 3 1 - 2K 2 + 2V 2 - W 2 K 2 - w 2 2 (w 2 - k) 2 

A simpler (the ordinary trace) prescription for the NBI Lagrangian consists in summation over color indices in the 
field invariants F^ v F a ^ v , F^ u F a ^ u in the square root form of the Lagrangian (|58|) . This gives 



L = -Aa'Vl - ZK 2 + 3V 2 - 9K 2 V 2 . (61) 



A. NBI cosmology 



Homogeneous and isotropic NBI cosmology with an ordinary trace Lagrangian turns out to be completely solvable 
by separation of variables [26| • It leads to an interesting equation of state: 

e(e c -e) 
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where e c — f3 /An is the critical energy density, corresponding to vanishing pressure. For larger energies the pressure 
becomes negative, its limiting value being 

p = -e/3. (63) 

This is the equation of state of an ensemble of non-interacting isotropically distributed straight Nambu-Goto strings 
(which indicates on the stringy origin of the NBI Lagrangian). In the low-energy limit the YM equation of state 
p = e/3 is recovered. Thus, the NBI FRW cosmology smoothly interpolates between the string gas cosmology and the 
hot Universe. The energy density is 

/ / a 4 + 3(w 2 -k) 2 \ 

e=ec (y « 4 -3^ (64) 

From the YM (NBI) equation one obtains the following evolution equation for the energy density: 

. = _ 2 a e(e + 2e c ) ^ 

a £ + £ c 

which can be integrated to give 

a 4 (e + 2e c )e = const. (66) 

From this relation one can see that the behavior of the NBI field interpolates between two patterns: 1) for large energy 
densities (e >• e c ) the energy density scales as e ~ a~ 2 ; 2) for small densities e«e c one has a radiation law e ~ a~ 4 . 
Remarkably, the equation for the scale factor a can be decoupled (g = (3G): 

2ga(a 2 + k) 



2ga 2 + 3(ct 2 + k) 
and admits the first integral 

3(a 2 + k) 2 + 4ga 2 (a 2 + k) =C, (68) 
which allows to draw phase portraits for different k 

• Closed. The only singular point is a = 0, b = 0, (b = a) which is a center with the eigenvalues ±i\/6g. Solutions 
evolve from left to right in the upper half-plane as time changes from — oo to oo, and from right to left in the 
lower half-plane. All solutions are of an oscillating type: they start at the singularity (a = 0) and after a stage of 
expansion shrink to another singularity. The global qualitative behavior of solutions does not differ substantially 
from that in the conformally invariant YM field model, except near the singularity: 

a ( t ) = b t-^t 3 + O(t 5 ), (69) 

where 6 is a free parameter. Absence of the quadratic term means that the Universe starts with zero acceleration 
in accord with the equation of state p w —e/3 at high densities. 

• Spatially flat. There is a singular line b — each point of which represents a solution for an empty space 
(Minkowski spacetime). This set is degenerate, and there are no solutions that reach this curve for finite values 
of a. All solutions in the upper half-plane after initial singularity expand infinitely. A remarkable fact is that for 
this case one can write an exact solution for a in an implicit form: 

Vff (*- *o) = V3(n - arctanfT 1 +tt/2) , (70) 

where = y/2a/ \/ \/ a 4 + ~C — a 2 . The metric singularity is reached at t = to. 

• Open. Physical domain is & < — 1, 6 > 1. There is a center at a = 0, b = with the eigenvalues ±iy/6g, but it 
lies outside the boundary of the physical region. Other singular points are (a = 0, b = ±1). These points are 
degenerate and cannot be reached from any point lying in the physically allowed domain of the phase plane. The 
only solutions which start from them are the separatrices b = ±1 that represent (part of) the flat Minkowski 
spacetime in special coordinates. One can easily see that all solutions in the upper part of the physical domain 
a > 1 start from the singularity and then move to a — > oo, a — > 1. 
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B. NBI on the brane 

Replacement of the standard YM Lagrangian by the Born-Infeld one breaks conformal symmetry, providing 
deviation from the hot equation of state and creating negative pressure. Surprisingly enough, putting the same NBI 
theory into the RS2 framework gives rise to an exact restoration of the conformal symmetry by the brane non-linear 
corrections j27|. Choosing the ordinary trace action 



^- det(<fc„+F MV //3) d A x -k 2 (R 5 + 2A 5 )V^d 5 x, (71) 



S = ATr 

where the brane tension A plays a role of the BI critical energy density, one obtains the constraint equation 



- = -A+-(A + £ ) 2 + ---, (72) 
a J 6 36 ar a z 

where £ is the integration constant corresponding to the bulk Weyl tensor projection ("dark radiation") and, as usual, 

A 4 = i« 2 (A + i. 2 A 2 ), G (4) = ^. (73) 
The energy density in this model scales as 

s = A (V(l + C/a 4 - l) , (74) 

where C is the integration constant. Surprisingly, the constraint equation comes back to that of the YM conformally 
symmetric cosmology 



2 

a ' 



~~ 3~~ A4 + ^"^' (75) 

where the constant C — £ + k 4 A 2 C/36 includes contributions from both the "dark radiation" and the YM energy 
density. 



V. CONFORMAL SYMMETRY BREAKING AND DARK ENERGY 

Conformal symmetry breaking in NBI theory demonstrates the occurrence of the negative pressure, but its extremal 
value p = — e/3 is still insufficient for DE. Meanwhile, a stronger violation of conformal symmetry may provide an 
equation of state with e ~ — 1 . Such violation can be of different nature: 

• Quantum corrections, 

• Non-minimal coupling to gravity, 

• Dilaton and other coupled scalar fields including Higgs, 

• String theory corrections. 

Here we just explore some model Lagrangians to see the necessary conditions for DE. Assuming the Lagrangian to be 
an arbitrary function L(J-,Q) of invariants 

T = -F% v F a » v /2 and Q = -F£ v F a » v /4, (76) 
one finds for the pressure and the energy density (conformal time): 

e = L + (&^w 2 + 3^w(k - w 2 )^) a-\ (77) 
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For a simple estimate consider the power-low dependence: 

L~F{F/ij?y-\ (78) 
where fi has the dimension of mass. Then in the electric (kinetic) dominance regime one obtains 

w =-=4p^ty < ra > 

e 6{2u — 1) 

For certain v this quantity may be arbitrarily close to W = — 1 or even less. An electric phantom regime is thus 
possible. 

In the magnetic (potential) dominance regime one obtains 

W = 4i//3 - 1. (80) 

The value W = — 1 can not be reached, but an admissible DE regime is also possible. These regimes are transient 
since during the evolution the electric part transforms to the magnetic and vice versa. 

Another plausible form of the lagrangian (suggested by quantum corrections) is logarithmic @ : 

L ~ Tln(T/n 2 ). (81) 
e = 3 (T[ln( J"//! 2 ) + 2] + V ln(J/> 2 )) a~ 4 , (82) 



Then the energy density is 



and the equation of state is 

w = p - = ( t + v)H^/» 2 ) + Wv~t) 

e 3(T + l/)ln(J-//! 2 ) + 6T ' [ ' 

where T — w 2 , V = (k — w 2 ) 2 . It is easy to see that W ~ —1 for ln(.F//i 2 ) ~ — 1. In this case the DE regime is 
transient. The phantom regime is also possible. 

Thus, the DE conditions for the homogeneous and isotropic YM field can arise indeed as a result of sufficiently 
strong breaking of the conformal symmetry. 



VI. MISCELLANIES AND OUTLOOK 



Though not directly related, we would like to mention here an interesting attempt to derive dark energy scale 
cosmological constant from gluon vacuum of QCD by Klinkhamer and Volovik [30j. The issues related to chaotic 
behavior of homogeneous non- isotropic YM configurations were discussed in [3l| . Possible applications of YM fields 
to cosmology up to partial identification of CMB with "cold" YM matter as discussed in the Sec. 2 were recently 
proposed by Tipler [32j . We want, however, to conclude with more conservative viewpoint that further work is needed 
to understand whether YM fields can be relevant in realistic cosmology indeed. First of all this includes deeper 
investigation of the coupled YMH dynamics and a thorough analysis of cosmological perturbations in various EYM 
scenarios. 
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Phc. 1: The amplitudes h(t)/a, f(t) of Higgs and YM fields and the state parameter w(t) = p/e for a typical solution with 
long quasistationary state. 
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Phc. 2: Phase trajectory for the scale factor in vicinity of the stationary state. 
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Phc. 3: Phase tracjectory of the Higgs field demonstrating the jump of the minimum of the scalar field potential from hmm — 
(equals to zero when j — 0) to a with the growing scale factor. 
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Phc. 4: The field amplitudes and the scale factor for two subsequent cycles with different parameters. 
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Phc 5: A typical solution for small a (on this plot a = 10 ). Initially the scalar field amplitude is large enough to ensure 
cosmic acceleration. With decreasing amplitude an accelerated expansion changes to contraction. 
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Phc. 6: The number of e- folds during inflation in the case of simple slow roll inflation, 7V_, and hybrid inflation, JV+, as functions 
of initial value of the scalar field, ho. 
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Phc. 7: The dependence of the number N of e- folds on the initial value of the scale factor oo, as compared with the critical 
value a c ~ \^6/(a/3). When ao reaches a c , the hybrid inflation reduces to the simple slow roll inflation. The peaks in the left 
part of the plot describe the divergence of N due to an infinite inflation, when the Higgs field is driven exactly to the top of 
the Higgs potential. 



